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Abstract 

Using Gelfand-Levitan theory on a half line, we construct a vessel for the class of potentials, 
whose spectral functions satisfy a certain regularity assumption. When the singular part of the 
spectral measure is absent, we construct a canonical model of the vessel. Finally, evolving the 
constructed vessel, we solve the Korteweg de Vries equation on the half line, coinciding with the 
given potential for t = 0. It is shown that the initial value for a; = is prescribed by this construction, 
but can be perturbed using an "orthogonal" to the problem measure. 

The results, presented in this work 1. include formulas for the ingredients of the Gelfand- 
Levitan equation, 2. are shown to be general in the sense that NLS, Canonical systems and many 
more equations can be solved using theory of vessels, analogously to Zacharov-Shabath scheme, 
3. present a generalized inverse scattering theory on a line for potentials with singularities using 
pre-vessels, 4. present the tau function and its role. 
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1 Introduction 



The Korteweg-de Vries (KdV) is the following nonlinear evolutionary Partial Differential Equation 
(PDE) for a function of two real variables q(x,t): 

qt = -^qq:, + ^qx^x, (1) 

where qt,qx denote the partial derivatives. The mathematical theory behind the KdV equation 
is rich and interesting, and, in the broad sense, is a topic of active mathematical research. The 
equation is named after Diederik Korteweg and Gustav de Vries who studied it in [DJK95| . Usually, 
one considers the initial value problem, which is defined as follows: find a solution q{x,t) of ([T]), 
which additionally satisfies: 

• q{x, 0) = q{x){x £ R), for a given function q{x) and is called KdV on line, or 

• q{x,0) = q{x){0 < x < oo), q{0,t) = h{t){0 < t < oo), called KdV in the first quadrant. 

A more complicated problem consists of specifying the initial value on a closed loop. A standard 
technique to solve KdV involves a more elementary equation, called Sturm Liouville (SL) differential 
equation: 

- -^y(^) + q{x)y{x) = Xy{x), (2) 

where A G C is called the spectral parameter and the coefficient q{x) is called the potential. In order 
to solve IT]) using ^ one transforms [CG67| the potential q(x) appearing in to its "scattering 
data" . Then one evolves with t the scattering data using some simple differential equations. Finally, 
transforming back the evolved scattering data we obtain a new potential q{x, t) of two variables, 
which solves IT} and satisfies 5(2;, 0) = q{x). In other words, in this manner we solve the initial 
value problem for the equation |l]). So, in order to solve |T]) one have to find "scattering data" for 
a given potential. 

On the half line the question of characterizing of scattering data (or more precisely spectral 
function) for a given potential was completely solved for a continuously differentiable potential by 
Gelfand-Levitan theory [IMG51) . but not always it can be used to solve the KdV equation |T]). In 
the case the spectral function is "differentiable" (var[p{X\ < cxa) it is possible to use it and this fact 
is reflected in many known solutions corresponding to this case: 

1. Soliton solutions correspond to dp{\) to be a finite sum of point mass measures (discrete 
measure), 

2. Krichever solution |Kri77| . where dp{X) = f{X)d\ for A G F for some algebraic curve F, 

3. Fadeyev inverse scattering theory [Fad74| where the dp(\) is supported on the positive real 
line and has a finite number of point-mass measures on the negative real line, 

4. Periodic potentials [Mag66| correspond to discrete spectral measures with accumulation point 
at infinity, 

and many more. In all these examples, the potential q{x) exists on the whole line R. Gelfand-Levitan 
theory |IMG51| , on the other hand, was developed primary to answer the question of reconstructing 
the potential q{x) from its spectral function p(A) on the half line using some works of Marchenko 
and H. Weyl. 

An analogous scheme for solution of the KdV equation |l| was developed by the author using a 
theory of vessels |Melc| . In this theory the scattering data is encoded (and easily translated) to a 
matrix- valued 2x2 function ^(A), which is analytic at infinity, with value I there and is symmetric 
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(see Theorem 12. 2p . As an initial step so called "regular" vessels were used, which involved only 
bounded operators. It is known |BGR90| that such functions possess realizations as follows 



1 
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id, 



where for an auxiliary Hilbert space H the operators act as follows: Xo,A ; H — )■ "H, 
Bq : ~^ H and satisfy the following conditions 



Evolving the operators with respect to x, one uniquely constructs 'Melal a vessel 5J, whose transfer 
function S{\,x) is realized as follows (|24|l : 



and coincides with the given S{\) for x = xq (see Definition 12.21 for details). The operator Xq 
is assumed to be invertible and A is assumed to be a generator of a Co-group. Notice that the 
singularities of 5(A) are determined by the spectrum of A only. Moreover, the same singularities 
with respect to A appear in S{X,x) and are a;-independent. The new function S{\,x) serves as a 
Backlund transformation (see Theorem 1 2. 5 P between trivial SL equation (having the zero potential) 

cP 

and a more complicated one, having potential 2—— ln(det(X(3;))). Evolving next the operators with 

dx^ 

d^ 

respect to t as in Definition 15 . 1 1 we will obtain a KdV vessel, for which q{x, t) — 2^^-^ ln(det(X(x, t))) 

satisfies As we can see this scheme is essentially the same as the Zacharov-Shabbath's one 

[ZABS74 ], but has no assumption on the (rapid) decrease of the potential at infinity. Moreover, one 
can use vessels to construct solutions for other completely integrable PDEs, for example NLS was 
done in Molb . 

In this work we present a generalization of regular vessels, involving an unbounded operator A. 
The main problem in this approach is that many formulas, present in a regular vessel, have to be 
rewritten, either including the domain of A, or using its resolvent -R(A) = (AJ — A)~^ . We have 
chosen the second approach, so for example, the Lyapunov equation 



If all the operators are bounded these two equations are equivalent. All the regularity assumptions 
and formulas are rewritten so that Gelfand-Levitan theory fits into the setting of vessels in a par- 
ticular case of the spectral measure. More precisely, we construct a vessel in the case when a new 
function if (A), constructed from the spectral measure p{X) 



AXo + Xo^l* + BoaiBo = 0, XS = Xq. 



S{X,x) = 1 - B' {x)X{x)~\\I - Ay^B{x)ai 



AXo + XoA* + Bo(JiBo = 



becomes in the new setting (I12|) 



XoR*{-\) + i?(A)Xo - R{X)BoaiBoR''{-X) = 0, A ^ spec(^). 




A > 0, 
A < 0, 



is a measure and has the property that (f8^ holds: 




4-00 



— oo 
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The function f{x, y) is one of the ingridients of the famous Gelfand-Levitan equation HlOf) . We show 
in this case how to construct the corresponding KdV vessel, and as a result a solution of the KdV 
equation ([T]) coinciding with the given continuous q{x) (whose spectral function creates /(a;, y) via 
(fST ). The value of the vessel for x = 0, t > turns to be prescribed and we present a formula for it 

(gHli 

q{Q,t) = 2 I y/\sm{{-s/\ft)du[\). 

It is also possible to perturb this value, without changing the value for q{0,x), which is explained 
in Section [5l 

Notice that the solution of the KdV equation q{x,t) = 2— — - ln(det(X(a;, t))) has the property 

ax^ 

that it does not exist at the points where X(a;, t) is singular. At these points the potential q{x, t) 
has a singularity, of course. So, we propose to study potentials, which arise in such a manner using 
a notion of pre-vessel (see Definition 14. 1|) . 

At the last Section [6] we show some models of vessels with an unbounded A and a realization of 
the Fadeyev theory. The canonical systems and Non Linear Schrodinger equations are presented as 
particular cases of the theory of vessels as a final example. 



2 Background 

2.1 Summary of the Gelfand Levitan theory 

A fundamental paper of Israel Gelfand and Boris Levitan |IMG51) serves as a background and a 
reference for many research works on SL equation on a half line, . Let q{x){Q < x < oo) be a 
continuous on each finite interval function. Consider the equation ([2} with boundary condition 

2/(0) = 1, y'{0)^h (3) 

and let (/>(a;, A) be its solution. There exists a representation of (j}{x, A) in the following form 



(j}{x,X) = cosiVXx) + J K{x,t) cos{y/x t)dt. (4) 



And conversely, there is a representation 

X 

cos{V\x) = <j){x,X) - J Ki{x,t)<j){t,X)dt. (5) 



From the work of Marchenko |Mc77| using Riemannian transformations, it follows that for contin- 
uously differentiable q{x) and h = 

^K{x,t) - q{x)K{x,t) = ^K{x,t), 

|i^(x,t)|..o, £i^(x,x) = /;^dx. 

The second kernel satisfies 

^K^{x,t) = ^K,{x,t) - q{t)K{x,t), 
|i^(x,t)l,.o, ^K(x,x) = !;^-^dx. 
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These PDEs with initial value problems are known to have unique solutions. 

There exists at least one non decreasing function p(A)(— oo < A < +00), such that for each 
f{x) £ L2(0,oo) functions 



E,r{X) = J f{x)4>{x,\)dx 


converge in norm with respect to p(A) to a function -B(A), i.e. 

+00 

lim / [E{\)-E„{\)]dp{\)^Q. 

n— >-oo J 

— OCi 

In this case the Parseval equality holds: 

+00 

E\\)dp{\) = j f\x)dx. 

—00 

The function p(A) is called the spectral function of ([2]) with initial conditions Q. It is proved in 
Gelfand-Levitan paper that p(A) satisfies the conditions of the next Theorem. 

Theorem 2.1 ( [IMG51] ). Let p(A) be a monotone function. Let 



and suppose that to (A) satisfies 
1° for all X > there exists integral 



p(A) - -VA, a > 0, 
p(A), A<0 



'■duj{X), 



2° the function 



possesses the fourth continuous derivative. 

Lf the function p(A) (or more precisely io{X)) satisfies conditions 1°,2°, then there exists a con- 
tinuous function q(x) and h, such that p(A) is the spectral function of ([2]) with initial condition 
Q. Conversely, if q{x) is continuously differentiable, the corresponding to it spectral function p(A) 
satisfies conditions 1°,2°. 

The smoothness of q{x) is determined by uj{X) as follows; if a{x) has derivatives of degree n + 4, 
then q{x) possesses n derivatives. Conversely, if q{x) has n + 1 derivatives, then a{x) has n + 4 

. — 

continuous derivatives. Notice also that J eV'^'^da;(A) defines an analytic function on (0,oo), 

— 00 

because its derivatives exists for all n: 



u 

/ iVwr 



"eV l^l-da;(A) < / eV l^l(-+")da;(A) < 00 
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following from \/|A[ < ev The existence of the derivatives of a{x) is directly related to the 
convergences of the integrals for all x G (0, oo) for consecutive values of n 

oo oo 

x)du}{\). 



One direction of Theorem 12.11 is classical and very well studied. The second direction of recon- 
struction of q{x) and h for a given p(A) involves the function 

F{x,y) = |°£MVM£in(VMrf,(;,), (7) 

— 00 

which turns to possess the second mixed derivative 

The authors notice that when lo{X) behaves "nicely" at infinity (if for example var[u}{X)\ < 00) 
fix, y) can be represented as follows 

/(a;,y) = ^^= J cos{V\x)cos{v^y)dLj{\). (8') 

— 00 

The potential q{x) turns to be 

<lix) = 2'-^, (9) 
where K{x,y), appearing also at Q satisfies the famous Gelfand-Levitan equation 

X 

fix, y) + Kix, y)+ f Kix, t)fit, y)dt = 0. (10) 



2.2 Vessels 

In a series of papers of the author [Mel09l iMelal IMelcl iMelb] and coUaborators [AM091EMT2] there 
was developed a similar approach to the inverse scattering of many LDEs (like SL, NLS, Canonical 
systems) using a theory of vessels. In the most general setting in the theory of vessels, the spectral 
function p(A) is translated into a complex-valued matrix function, belonging to a special class of 
functions, parallel to the "scattering matrix" s(A) appearing in the work on Inverse scattering of 
Faddeyev [Fad74) : 

Definition 2.1. Let ai = al be an invertible pxp matrix. ClassUiai) consist ofpxp matrix-valued 
functions S'(A) of the complex variable A, possessing the following representation: 

Six) = I - BoXa\XI ~ Ay^Boai (11) 

where for an auxiliary Hilbert space H there are defined operators Bq : C H, A,Xo : H H. 
A general matrix-function S'(A), representable in such a form is called realized. Moreover, the 
operators are subject to the following assumptions: 
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1. the operator A has a dense domain D{A) with a symmetric spectrum with respect to the 
imaginary axis. A is the generator of a Co semi-group on T-L. Denote the resolvent of A 
as R{\) = {XI - A)~^ , 

2. the operator Bo satisfies R{\)Boe G H for all X ^ spec(^), e £ C. 

3. the operator Xo is bounded, self adjoint, invertible, 
4- the Lyapunov equation holds for all X ^ spec(yl); 

XoR'i-X) + R{X)Xo - R{X)BoaiB;R*{~X) = 0. (12) 

In this case S{X) is also called symmetric. 

We call the last equation Lyapunov, because in the case all the operators are bounded, we obtain 
a usual Lyapunov equation AX + XA* + BociBq ~ 0. It is easy to show using the realization pip 
and the Lyapunov equation only that the symmetry condition S* {—X)aiS{X) — ai holds for all 
A ^ spec (A): 

Theorem 2.2. Suppose that S(X) G then S* ("A)criS'(A) = cti for all X ^ spec(yl). 

Proof: Let us plug first the realization formula pil) for ^(A) into the symmetry condition: 

5'*(-A)aiS(A) = (J + aiB*o{XI + A'^^Xg^ Bo)ai{I - B^X^^XI - Ay^Boai) 

= (71 + aiB^iXI + A*)-^X^^Boai - aiB^H^^XI ~ A)-^Boai~ (13) 
-aiB^iXI + A'y^Xg^BoaiB^X^^XI - A^^Boai 

To the last term we apply the Lyapunov equation, rewriting it: 

aiBo*(A/ + A*)-^Xo ^BocriBSXo i(A/ - A)-^Boai = 

= aiBoiXI + A'y^Xo^iXI - A){XI - A)-^BoaiBo{XI + A*)-^{XI + ^*)X,7 ^(A/ - A)-^Bocri 
= aiBolxi + ^*)-^Xo ^(A/ - A)[Xo(A/ + A*)-^ - {XI - Ay^Xo\{XI + A*)Xq^{XI - Ay^Booi 
= criBo*(A/ + 4*)-'Xo ifiocri - aiBSXo ^(AJ - A)-^Bqg^ 

If we plug this expression back into (|13|) we obtain the symmetry condition. □ 
When S(A) is just analytic at infinity (hence A must be bounded), there is a very well known 
theory of realizations developed in [BGRQOj . For analytic at infinity and symmetric, i.e. satisfying 
S*(— A)criS'(A) = CTi, functions there exists a good realization theory using Krein spaces {H is a 
Krein space), developed in [ADp l. Such a realization is then translated into a function in 14 (ai). In 
order to acquit the name of "scattering data" for S{X) we show how to construct the "potential" 
(analogously to SL case), which we call 7* (a;). For this we will need to fix two additional p x p 
matrices = (T2 and 7 = —7*: 

1. Suppose that there exists a solution -8(2:) satisfying B{0) = Bo and the equation 

0= -^(R(X)B(x)aie) + AR(X)B(x)(j2e + R(X)B(x)'ye, VA ^ spec(A), e G C''. (14) 
ax 

This equation is solvable because the coefficients cti, 0-2,7 constant and A is a generator 
of a Co semigroup. We will solve it in an exact manner for the SL case. For this equation to 
hold we also demand that the following regularity assumptions hold: 

VA ^ spec(A) : R{X)B{x)a2e' C D{A), (15) 
VA ^ spec(^) : R{X)B(x)'yC' C H. (16) 

^At the paper jADI a similar result is proved for functions symmetric with respect to the unit circle, but it can be 
translated using Galley transform into 5'*(— A)ctiS'(A) = cti and was done in IMelal IAM12| 
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2. Solve for X{x), satisfying X(0) = Xo and: 

-^X{x) ^ B{x)a2B*{x). (17) 

Notice that B{x)a2C^ C H follows from (fTS)) . 

3. Define "f*{x) on I by 

-y^ix) = 7 + a2B*(x)X~^{x)B{x)ai ~ aiB*{x)X~^{x)B{x)a2. (18) 

The main reason, why we call S{X) as the "scattering data" is the fact that 7* (a;) (generalized po- 
tential) is uniquely determined from S{X) by this construction. With some additional requirements, 
there is also a unique S{X) for each 7* (a:). Still, the question of reconstructing 5(A) from the given 
7«(a;) is the fundamental problem of the theory of vessels, solution of which we are going to present 
here for the Sturm Liouville case under a certain regularity assumption (Theorem I3.4|l . Moreover, 
the function 

S{X,x) = I ~ B*{x)X~^{x)(XI - A)''^B{x)ai (19) 

acts as a Backlund transformation, because (see Theorem I2.5|l multiplication by it maps solutions 
m(A, x) of the input LDE with the spectral parameter A: 

d 

-ai—u{X,x) + {a2X + 'y)u{X,x)^0 (20) 

to solutions y{X, x) = ^(A, x)u{X, x) of the output LDE with the same spectral parameter: 

- ai^y(A,x) + {a2X + 'y4x))yiX,x) = 0. (21) 

Notice also that we relate in this manner solutions of an LDE with constant coefficients to solutions 
of an LDE with variable coefficients analogously to the scattering theory. 
Such functions are called transfer functions of vessels, defined in the next 

Definition 2.2. The collection of operators and spaces 

2J = (A, B(x),X(x); ai, aa, 7, 7,(a;);1H, I), (22) 

is called a vessel, iflL(x),A -.T-i-^'H, B{x) : C — >■ H are differentiable linear operators, subject to 
the following conditions: 

1. the operator A has a dense domain D{A) with a symmetric spectrum with respect to the imag- 
inary axis. A is the generator of a Co semi-group on H, 

2. B(x) satisfies regularity assumptions (|15|) . (|16|l . and the equation (|14[) . 

3. X(a;) is bounded, self adjoint, invertible on I and satisfies H17|) . 
4- the Lyapunov equation holds for all x £ 1, X ^ spec(yl); 

X(2:)i?*(-A) + R{X)Xix) ~ R(X)B{x)aiB(x)*R''{~X) = 0. (23) 

5. 7.(2;) satisfies (|18|) . 

The class of transfer functions is defined as follows 

Definition 2.3. Class X(a"i) consist of p x p matrix-valued (transfer) functions S{X,x) of the 
complex variable X and x £l for an interval I — [a, b] possessing the following representation: 

S{X,x) ^ I ~ B*{x)X'^{x){XI - Ay'^B{x)ai (24) 

where for an auxiliary Hilbert space T-L, the operators B{x) : C — ?> T-L, 'K(x),A : H ^ H constitute a 
vessel for some (72,7* (a;) on I. 



8 



we want to prove the Backlund transformation theorem, presented earUer. Before we do that 
few techincal results are needed: 

Lemma 2.3 (permanence of the Lyapunov equation). Suppose that B{x) satisfies (|14p and 
X(a;) satisfies (|17p . then if the Lyapunov equation ((23} holds for a fixed xo, then it holds for all x. 
//X(a;o) = X*(a;o) then X(a;) is self-adjoint for all x. 
Proof: Let us differentiate Lyapunov equation (1231) : 

^[-X(a;)(A7 + A'y^ + {XI - A)-^X{x) + {XI - A)'^ B{x)aiB'' {x){XI + A*)-''] = 
= -B{x)a2B*{x){XI + A*)-'' + {XI - A)-^B{x)a2B*{x) + 

+ ^[(A/ - A)-^B{x)aiB''{x){XI + A")-^] = 
= -B{x)a2B'-{x){XI + A*)-'^ + {XI - A)-'^B{x)a2B'-{x) + 

-A{XI - A)-'^B{x)a2B*{x){XI + A*)-'^ - {XI - Ay^ B{x)iB* {x){XI + A*)-'^- 

-{XI - A)-^B{x)a2B'-{x){XI + A*)-^A'- - {XI - A)-^ B{x)/'B* {x){XI + A')"^ 

The terms involving 7 are canceled, because 7 + 7* = 0, by the assumption on it. So this becomes, 
after adding and subtracting XI next to A and to A*: 

= -B{x)a2B*{x){XI + A'y + {XI - A)-^B{x)a2B* {x)+ 
-{A -XI + XI){XI - A)-^B{x)a2B*{x){XI + A*)''^- 

-{XI - A)-'^B{x)a2B*{x){XI + A*)-'^{A* + XI- XI) = 
= -B{x)a2B''{x){XI + A*)-'^ + {XI - A)-'^B{x)a2B*{x) + 

+B{x)a2B'-{x){XI + A*)-^ - X(AJ - A)-^ B{x)a2B* {x){XI + A*)-^ 

-{XI - A)-^B{x)a2B''{x) + {XI - A)-^ B{x)(j2B'' {x){XI + A*)-^^ = 

= 0. 

So, it is a constant and the result follows. 

Since the derivative (I17p — X(a;) = B{x)a2{x)B*{x) is self-adjoint, X(x) will be self-adjoint, 
dx 

once X(a;o) is. □ 
Lemma 2.4. Suppose that B{x) satisfies p4|) and X(a;) satisfies then 

^[{XI + A*y^X-\x)B{x)ar] = 
dx 

- X{XI + A*y^X-^{x)B{x)a2 + X~^{x)B{x)a2 - {XI + 4*)"^X"^(a;)B(a;)7,(a;). (25) 
Remark: 1. Notice that if A is bounded it is equivalent to 

■^[X~^B{x)ai] = A''X~^B{x)a2 - B{x)-y,{x) (26) 



2. Substituting A with —A in (I25|) and taking adjoint, we will obtain that the following equality 
holds 

-^[a^B*{x)X-\x)R{X)] = 

Xa2B''{x)X~^{x)R{X) - a2B* {x)X-^ {x) - j4x)B* {x)X-^ {x)R{X). (27) 
Proof of Lemma [2T4l Let us differentiate {XI + A'y^X''^ B{x)ai using ([14}, pT} : 

^[{XI + A')-''X-'B{x)ai] = 

= {XI + yl*)-i-^[x-ilB(a;)ai + {XI + ^*)-iX-i(A/ - A)-^[{XI - A)-^B{x)ai] = 

dx dx 
= -{XI + A'')-^X-^B{x)a2B*{x)X~^B{x)ai + 

+(AJ 4- A*)-^X-\XI - A)[-A{XI - A)-^B{x)a2 - {XI - A)-^B{x)'y]. 
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Notice that AJ — A is a bounded operator for A ^ spec (A). Then rearranging 

= -(A/ + A'y^X-^ B{x)a2B'' {x)X~^ B{x)ai - (AJ + A*)"^X"^B(a;)7- 
-(A/ + ^*)-^X-^(A/ - A)A{XI - A)-'^B{x)a2. 

If we add and subtract AJ next to A at the last term, we obtain 

= -(A/ + ^*)-^X-^B(x)o-2B*(a;)X"^B(a;)cri - (AJ + A*)-'^X-^ B{x)'y- 
-(AJ + A*)-^X-^(AJ - A)(A - AJ + AJ)(AJ - Ay^B{x)a2 

= -(AJ + A*)-'^X-^B{x)[-y + a2B*{x)X-^B{x)ai]- 

-A(AJ + A*)-'^X-^B(x)a2 + (AJ + ^*)-^X-^(AJ - A)B{x)a2 



(28) 



A simple consequence of the Lyapunov equation (|23|) is as follows 

(AJ + yl*)"^X"^(AJ -A) = X-^{x) + (AJ + A'y^X'^ {x)B{x)aiB{x)X-^ (x) 
So if we plug this expression into the last term of ((28}, we obtain rearranging 

^[{\I + A*y^X-'B{x)ai] = 

ClX 

= -(AJ + A'')-'^X-^B{x)['y + a2B*{x)X~^B{x)ai]- 

-A(AJ + A*)-^X-^B{x)(j2 + [X-^{x) + (AJ + A'')-^X-^{x)B{x)aiB{x)X'\x)]B{x)a2 = 
= -(AJ + A*)-''X-^B{x)['y + <J2B* {x)X-^ B{x)cji - aiB{x)X-\x)B{x)a2]- 

-A(AJ + A*)-'^X-^B{x)a2 + X-^{x)B{x)a2 
= -(AJ + A'')-^X-^B{x)'y,{x)- 

-A(AJ + A'-)-^X-^B{x)a2 + X-^x)B{x)a2, 



using (|18|) at the final step. □ 
Now we have all the ingredients of the following Theorem. This theorem has its origins at the 
work of M. Liv§ic [LsOT] and was proved for bounded operators in [Melal iMelbl IAM12] . Now we 
present a generalization of these results for the case of unbounded operator A. 

Theorem 2.5 (Vessel=Backlund transformation). Let V be a vessel defined tn (1221) and satisfying 
all the conditions of Defimtion \2.2l Fix A ^ spec(yl) and let m(A, x) be a solution of the input LDE 
(|20p . Then the function y{X,x) = S{X,x)u{\,x) is differentiate and satisfies the output LDE \21\ . 



Proof: Let us fix A ^ spec(yl) and a solution u(\,x) of (|20|) . Then for y{X,x) = ^(A, x)u(A, a;) we 
calculate: 

ai-^y(X,x) = ai-!^[(I - B*(x)X-^x)(\I - A)-'^B(x)ai)u{X,x)] = 
ax d d 

= o-i— w(A,x) - ax-r\B*{x)X'^(x)(\I - Ay^B(x)aru(\,x)] 
dx dx 

= {a2\ + i)u(\,x)-Gi — [B*{x)X-^(x)] (XI - Ay^B(x)GM\x) 
dx 

-GiB*{x)X-^{x) — [(\I - A)-^B(x)]u{\x) 

-aiB*(x)X-^(x)(\I - A)-^Blx)ai4-u{\,x). 

dx 

Using ((57} and ((Ti)) it becomes 
ai^y{X,x) = {a2\ + "/)u{X,x) — 

--[Xa2B''{x)X-\x){XI - A)-^ - aiB'ixJX-^x) - 'y,{x)B* {x)X-\x){XI - A)-'^]B{x)aiu{X, x) + 
+aiB*{x)X-^{x)[A{XI - A)-^B{x)a2 + (AJ - A)-'^ B{x)^]u{X, x) 
-aiB*{x)X-\x)iXI - A)''^B{x){a2X + i)u{X,x) = 
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Let us combine the last two terms and rearrange the rest: 



= (a2A + 7)m(A, a;) — 

-\a2B*{x)X-'^{x){\I - A)-^B{x)aiu{X,x) - a2B*{x)X-^{x)B{x)alu{\,x)' 
'y4x)B'•{x)X-^{x){XI - A)-'^]B{x)aiu{X,x)+ 

+aiB''{x)X-^{x){A - A/)(AJ - A)-^ B(x)a2u{X, x) 
= (J2X(I - B''{x)X-^{x){XI - A)-'^B(x)ai)u(X,x) + 
+ [7 + cr2S*(a;)X"^(a;)_B(a;)cri - aiB*{x)X-'^{x)B{x)a2]u{X,x)~ 
-'y4x)B*(x)X-^(x)(XI - A)-'^B{x)aiu{X, x) 

Using p8[) and the definition of ^(A, a;) we obtain that it is 



= cr2A5'(A, x)u{X, x) + 7,(s)u(A, a;) — 
-'y,{x)B*{x)X-^{x){XI - A)-'^B{x)aiu{X,x) = 
= cr2A5'(A, a;)M(A, x) + 7,(a;)S'(A, x)u{X, x) = 
= (cr2A + 7,(a;))S(A, x)u{X, x) = 
= (cr2A + 7,(a;))y(A,2;). 

□ 

We have seen that there is a "simple" construction, called standard, of functions in I(cri), starting 
from a function in U{ai) by solving differential equations. Notice also that in the case Xq = / (the 
identity operator) we have a very well developed theory |MB58j of so called "characteristic functions" 
of the operator A, which has a finite dimensional real part: A + A* — ~B(t\B* . But in general, 
when Xo is not identity, it is only known that A has a spectrum, symmetric with respect to the 
imaginary axis. Existence of the vessel and its transfer function relies on the invertability of the 
operator X(a::). In order to investigate the existence of the inverse for X(a;) notice that from p7|l 



X(x) =Xo + / B{y)a2B' {y)dy 

it follows that 







So'X(a;) = J + Xo Y B{y)a2B' {y)dy. 



Since 0-2 has finite rank for dimf < 00, this expression is of the form I + T, for a trace class operator 
T and since Xn is an invertible operator, there exists a non trivial interval (of length at least ) 

iixo'ir 

on which X(a:) and t{x) are defined. Recall [IG69| that a function F{x) from (a, b) into the group 
G (the set of bounded invertible operators on H of the form I + T, for a trace-class operator T) is 
said to be differentiable if F{x) — I is dijferentiable as a map into the trace-class operators. In our 
case, 

-^(Xo 'X(3;)) = Xo '-^X(x) = Xo^B{x)a2B'-{x) 
ax ax 

exists in trace-class norm. This leads us to the following 

Definition 2.4. For a given vessel 2J ((22} the tau function t{x) is defined as 

T = det(Xo 'X(a:)). (29) 
The vessel exists at a neighborhood of every x for which r{x) 7^ 0. 
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Israel Gohberg and Mark Krein [IG691 formula 1.14 on p. 163] proved that if Xq ^X(a::) is a 
differentiable function into G, then t{x) = sp(X,7^X(x)) is a differentiable map into C* with 

= sp((Xo ^X(a;))-'^(x-ix(x))) = sp{X{x)'X-\x)) = 

= sv{B{x)a2B''{x)X~^{x)) = ^{023' {x)X'^{x)B{x)). (30) 
Differentiating this expression, we obtain that 

(-)' = — - {-f = ^t^{a2B''{x)X-\x)B{x)). 
T T T ax 

Using vessel conditions, since B{x), X~^{x) are differentiable bounded operators in the case AK~^{x)B{x) 
exists, or in the case it is canceled (SL case) we obtain that 

-^tr{a2B*{x)X-\x)B{x)) = 

= tT{a2a-^{-a2B*{x)A* - -y* B* {x))X-^{x)B{x))- 

-tr{a2B* {x)X-\x)B{x)a2B'' {x)X-^{x)B{x)) + tT{a2B'' {x)X-^{x){-AB{x)a2 - B{x)j)a^^) = 
= tr{a2a^^a2B*{x)[-A*X-^ - X-M]B(a;)) - tr([a2crj"S* + lajTVaJB* (x)X-^(a;)B(x))- 

-tT{a2B'{x)X-\x)B{x)a2B'{x)X-\x)B{x)) = 
= tr{a2a-\2B*{x)X-^B{x)aiB*{x)X-^B{x)) - tr([(72(7j:S* + 'y(7i^cT2]B* {x)X-^{x)B{x))~ 

-tr{a2B* {x)X-^{x)B{x)a2B* {x)X-^{x)B{x)). 

(31) 

In the special case of SL vessel parameters, we obtain that equations (|20p . (|2ip are equivalent to 
([2]). Let us explain it in more details. 

Definition 2.5. The Sturm Ltouville (SL) vessel parameters are defined as follows 





■ 


1 ■ 




1 


■ 




■ ■ 




1 













, 7 = 


i 



Suppose that we are given a SL vessel 2J, in other words, 03 (|22|) is defined for the SL vessel 
parameters. Denote a differentiable 2x2 matrix function ffo = B* (xy%r^ (x^Bix) = 



d 



Then from (O it follows that — = tx{o2B* (xy%r^ (x)B{x)) = a vs. the 1,1 entry of Hq. Using 

r 

simple calculations it follows that 

020-1^02 = 0, tT{o2Ho02Ho) = {ti{a2Ho))^ , o-20-^^-y* + fo^^a2 = 
and we obtain from (I31II that 



-i 

1 



= -tr( 



-i 

1 



Ho) = i(6* - b). 



Notice that the terms involving operator A are canceled. Moreover, we obtain that 

// / 

.T T 

h — h* n 

j,{x) = 7 + 02Ho(Ti — aiHo02 = 



b — b* a 
—a i 



r, r 

T 

4 

r 



Thus we obtain the following lemma (appearing already in [Mel09l Proposition 3.2 ]) 



^sp - stands for the trace in the infinite dimensional space. 
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Lemma 2.6. For SL vessel parameters, the following formula for "f*{x) holds 

7.(2;) = 



I 

r 



Analogously to |Mel09l Section 3.1.1], simple calculations show that denoting ?i(A, s) 
we shall obtain that the input compatibility condition pop is equivalent to 



ui(A, a:) 
U2(A, a;) 



The output y{\,x) 
equivalent to 



y2{\x) 



M2(A,a;) = -i-§^ui{\,x). 



S{X,x)u{X,x) satisfies the output equation (pT|) . which is 



-■§^yi{X,x) - 2— [ln(r(a;))]yi(A,x) = -i\yi(X,x), 



i/2(A,x) = -i[-§^yi(X,x) + —yi(X,x)]. 

Observing the first coordinates ui{X,x),yi{X,x) of the vector-functions u(X,x),y{X,x) we can see 
that multiplication by ^(A, x) maps solution of the trivial SL equation (i.e. q{x) = 0) to solutions 
of the more complicated one, defined by the potential 



q{x) 



-2^[ln(r(x))], 



(32) 



which can be considered as an analogue of the scattering theory. 



3 Construction of a vessel using Gelfand-Levitan theory 

Throughout of this section we suppose that /i = at the initial value problem ([Sjl . This will enable 
us to use the formulas @ and (O . Suppose that we are given a potential q{x) possessing a spectral 
function p{fJ,), for which f{x,y) satisfies the condition (f8^ : 



f{x,y) = 



dxdy 



= / cos 



i{y^x)cos{y^y)dui{fi). 



From here it follows (for x = y) that cos{y/JIx) is in the space 

n = Ll = {gill) I j j<7(M)|'rf^(M) < 

K 

and the functional [cos{\fXx)\* , defined by this function: 

+00 

[cos{^x)]* g{ii) = I cos{y^ x)g{ii)duj{fi) 

— 00 

is well defined on "H. Moreover, by formula Q we also obtain that since the kernel K{x, t) is 
continuous, the functional [(j){x, A)]* is well defined. 
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We are going to present a construction of a vessel 2J realizing the potential q{x). Let us define 
the operator ^ = iA as the multiplication operator on iX. It is usually an unbounded operator with 



the obvious domain. Define nexll 

B{x) = [ cos(y^a;) i^/Jlsin{^/JIx) ] (33) 

Then for A spec(j4), (XI—A)^^B(x) = — ^ — f cosi^/JIx) i^/JIsin(^/JIx) 1 . It is almost readable 
A — 

that the conditions (|I5|I . (|16|l are fulfilled and the equation (|14|l hold: 

^(AJ - A)-^B{x)ai = (A/ - Ay^ [ iAcos(VAt) -V\sm{VXx) ] = 

= -^(AJ - Ay^B{x)a2 - (A/ - Ay^Bix)^. 



As a result, we obtain that 



/(x,y)= [ 1 ]B'{x)B{y) 



[10] B*(a;)Xo'B(y) 



(34) 



Define next operator X(x) -.T-i^'H using an arbitrary g £%. Here we present how to calculate the 
resulting function '%.{x)g at point A: 

X X 

(X(x)fl)(A) = + j Biy)a2B*{y)dy]g = 3(A) + j cos(VAy) j cos{^y)gi),)dLu{fi) (35) 
a 

notice that this operator is of the form I + P for a positive, trace-class operator P, so it is invertible 
on [0, 00). From these formula we obtain an explicit representation for K(x,y): 

Lemma 3.1. The following equality holds: 

+00 

K{x,y) = - f (?!>(x,A)cos(\/Ay)dLj(A). (36) 



where (j){x,\) is the solution of ([2)) satisfying (f){Q,X) — 1, ^'(O, A) = (|3]). 
Proof: Plugging the expression Q for 0(a;, A) we obtain 

+ 00 +00 X 

J (l){x, X) cos{-\/Xy)doj(X) = J [cos{^/Xx) + J K{x,t) cos{\/^t)dt]cos{\/^y)duj{X) 

— 00 —00 

+ 00 +00 X 

= J cos(v'Aa;) cos(v'A i/)dLj(A) + J J K{x,t) cos{\/^t)dtcos{\/^y)duj{X) = 

— QC — 00 

X +00 

= f{x,y) + J K{x,t) J cos(\/At)dt cos(VA j/)da;(A) = 

-00 

^ f{x,y) + J Kix,t)f{t,y)dt^ 


= -K{x,y), 



by the Gelfand-Levitan equation (|10p . □ 



^for negative A we understand cos(\/A) as cosh(^/|A[) and similarly for the sinus 
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Corollary 3.2. The following formula holds 



'%.{x)(j){x,X) = co6{-\/\x) = B* [x) Q 

Proof: Immediate from (|35|l and the formulas (I36|l . Q. 

It turns out that we can find the inverse of X(a;) as follows. 

Definition 3.1. The operator Y{x) : "H — >■ "H is defined as follows 



(37) 
□ 



(Y{x)g){\) ^ g{\) ~ / 0(y,A) / M)<?(M)rf^(p)d2/ 



(38) 



Notice that if we change the variables in K{x,y) we will obtain a Kernel, which satisfies the 
defining initial value PDE for Ki{x,y). By the uniqueness of the solution of this problem it follows 
that K\{x,y) = K{y,x). This enables us to show that Y(a::) is the inverse of the operator X(a;). 
Similarly to the calculation of X(a;)(/>(x, A) we obtain, using ([5]) that 



Y(a;)cos(%/Ax) = Y{x)B'{x) 



= 0(k,A). 



(39) 



Since as we mentioned earlier [(/)(x. A)]* is a functional on T-L the following formula holds (Y*(a:) 
Y(a;) from the definition) 

[<^(:e,A)]* = [10 ]B{x)Yix). 
Moreover, using this notation, we also obtain that 

^[Y(x)] = -0(x,A)[0(a;,A)]*. 

Let us perform a simple calculation for the derivative of X(a::)Y(a;): 



d_ 

dx 



[X(a;)Y(2;)] 



-f [X(^)]Y(x)]+X(x)-f [Y(x)] 
ax , , ax 



= B{x) 

= B{x) 
= 0. 



1 ]B*{x)Y(x)~%{x)(t)(x,\)[(t){x,X)Y 



[<j,{x,X)Y -B{x) 



[<j}{x,\)]* = 



And since X(0)Y(0) — II = I is the identity operator, we obtain that ^{x) is onto, and since it is 
of the form I + C for a compact operator C, it is also invertible by a Fredholm theorem, with the 

inverse Y(a;). We could alternatively prove that — [Y(a;)X(x)] = to deduce the same result. 

Theorem 3.3. The operator Xix) is invertible for all x > and its inverse is the operator Y{x), 
defined in (|38|l . 



Let us take 



K'{x,y) = -[ 1 ]B*{x)X'\x)B{y) 



(40) 
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Then Gerfand-Levitan equation pOI) holds 
K'(x, y) + fix, y) + J K{x, t)f{t, y)dt = 



^K'{x,y) + f{x,y)^][ 1 ] B* {x)X-\x)B{t) 



1 ] B*{t)X^,'B{y) 



dt 



= K'{x,y) + f{x,y)-[ 1 ] B* (x)X'\x) J B{t)a2B* it)dtX-' B^y) 



= using vessel condition H14p — 

= K'{x,y) + f(x,y)-[ 1 ] B* (x)X-\x){X{x) ^ Xo)X^' B{y) 



= K'ix,y) + f{x,y)-[l ] B* (x)X^' B{y) +[l ] B* (x)X-Hx)B{y) 



= 0. 



By the uniqueness of the solution of (1101) we conclude that K'{x,y) = K{x,y) and the potential 
q{x), which is obtained from K{x,x) using ((9]) is 



- 2^ ^ ^-2±[1 ] B*ix)X'\x)Bix) 



dx 



dx 



= -2A,[ln(^(.))], 



(41) 

which coincides with ((32}. In other words, we obtain the following 

Theorem 3.4 (Realizing a potential on the half line). Let p{X) be the spectral function of the 
potential q{x) for which f{x,y) is defined usmg (f8^ . Define a vessel 

2J = {A, B{x),X{x); 0-1, aa, 7, lI,C^; [0, 00)), 

obtained by the vessel construction, applied to A — i\, Xo = /, Bo = [ 1 ] . Then this vessel 
realizes the given potential q{x). 

Proof: the differential equations (I14|) . p7|) are satisfied by the construction. The Lyapunov equation 
(|23[) for a; = is satisfies, since then AXq + XqA* = for Xq = / and A - anti selfadjoint by its 

1 



definition. Notice also that Bo ~ 







so that BqoxBq — and as a result (I23|l holds. By 



Lemma 12.31 it will hold for all x by the permanency of the Lyapunov equation. The operator X(a;) 
is invertible for all x £ [0,oo) with the inverse Y(a;), defined in psp by the arguments proceeding 
the Theorem. The potential of the vessel coincides with the given potential q{x) by (|41|l . □ 



4 Generalized scattering theory on a line 

A closer examination of the previous results, which can also be seen from the GL theory, implies 
that the operators B{x),X{x) are globally defined for all a; G R. Indeed the formula (|33p 



and dM 



B(x) = [ cos{y/Jix) i^/Jls\n{^/Jix) 



X 

{x) = I + j B{y)a2B*{y)dy 



define operators for all a; £ R. This implies the following definition 
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Definition 4.1. The collection of operators and spaces 

pre03= (A,B(a;),X(x);ai,<72,7;H,C''), (42) 

is called a pre-vessel, ifX{x),A : H H, B{x) : C — >■ H are dijferentiable linear operators, 
subject to the following conditions: 

1. the operator A has a dense domain D{A) with a symmetric spectrum with respect to the imag- 
inary axis. A is the generator of a Co semi-group on H, 

2. B{x) satisfies regularity assumptions (|15|) . (|16|l and the equation (|14|l for all x £9., 

3. X(x) is bounded, self adjoint for all x £ M. and satisfies (|17|l and there exists a point xo G R, 
in which X(a;o) is invertible, 

4- the Lyapunov equation (|23p holds for all x £9., A ^ spec(A). 

Tau function r(x) — det(X~^(a:o)X(a:)) is well defined now and exists on R. Moreover, if r(a') = 
t(6') = and t{x) 7^ for a' < x < b' , for arbitrary I = [a, b], where a' < a < b < b' we obtain a 
vessel on I. As a result, we can make the following definition 

Definition 4.2. The generalized potential of a pre-vessel pteQJ is the function 7,(2;), defined by 
psp for all points where t(x) 7^ 0. 

As we can see for each pre-vessel pre2J there is a unique generalized potential, realized by it. So, 
a very important question arises, is whether each generalized potential 7,(2:) arises in this form. It 
is also important to classify the classes of such generalized potentials and to relate them to classes 
of pre- vessels. 

Problem 4.3 (Generalized scattering theory). Characterize the class Q'P of generalized potentials, 
arising from pre-vessels. 

Problem 4.4 (Generalized inverse scattering). Given a generalized potential q{x) G Q'P, find a 
pre-vessel realizing it. 

For example, using Gelfand-Levitan theory we obtain: 

Theorem 4.1. Let q{x)(0 < x < 00) be a potential for which the function f{x, y) possesses the form 
(|8^ . Then there exists a pre-vessel, whose generalized potential coincides with q{x) on the positive 
real line. 



5 Solution of the KdV equation in the first quadrant 

Let as consider the initial value problem of ((T)). Suppose that for a given q{x) we have constructed 
a vessel 03 as in Theorem 13.41 Let us evolve the operators of the vessel 03, using the following 
conditions: 

^^B{x,t) = iA^B{x,t), B{x,0) = B{x) (43) 

^X^iABasB* -iBazB'A* +iB'yB*, X(a;, 0) = X(x). (44) 

For these equations to hold we have to add some regularity assumptions: 

B{x,t)a2C' C D{A), B{x,t)'yC'' CH. (45) 

Notice that since the operator A = iX is "diagonal" in our setting the equation (|43p is a sort of wave 
equation, whose solution is very easy to create: 

B{x,t) = [ cos{V\x - (Vxft) iV\sm{V\x - (Vxft) ] (46) 

Here x,t>0 and cos(\/A) = cosh(A/jA[) if A < as usual in Gelfand-Levitan notations. 
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Theorem 5.1. Let B{x,t) be a solution of the partial differential equations H14[) . (|43p . Then the 
formula for 'K{x, t) is as follows 

X{x,t)=Xo+ [ B(y,t)cj2B*{x,t)dy+ 
Jo 

+ f [iAB{0,s)a2B*{0,s)-iB{0,s)a2B''{0,s)A'' +iB{0,sy/B*{0,s)]ds (47) 



Proof: The first equation (|17|l is immediate, since tiie last term depends on t only. In order to 
prove (|44|l . we can calculate using (I14|l and (|43p that 

^ t)a2B* {y, t)] = ^ [iAB{y, t)a2B* {y, t) - iB{y, t)a2B* {y, t)A* + iB{y, t)'yB* {y, t)] 



and the result follows by applying J dy to both sides. □ 



Definition 5.1. The KdV evolutionary vessels is defined as follows 

QJ = {A,B{x,t),X(x,t);ai,cr2,'y,'y*{x,t); L'i,C^;l), 

where A = ijj. and B{x,t),X{x,t) are defined by (|46l) . (1471) and satisfy the regularity assumptions 
(HI). 



Notice that for each fixed t, we obtain a vessel, for which all operators depend on t and realize a 
potential q{x,t). For the KdV equation ^ we have to require the existence of the third derivative 
with respect to x. In this case, it turns out that the regularity assumptions H45[) imply that the 
expressions 



1 ] B'ix, 



\x,t)A'B{x,t) 



i = 0,1 



are meaningful. This enables to use the theory of evolutionary vessels, or more precisely [Melcl 
Theorem 10]: 

Theorem 5.2 (Vessel^ Global Solution in First Quadrant). Let ^Kdv be a KdV vessel, realizing a 
given potential q(x), defined on the half line with 4 continuous derivatives. Then the potential q{x, t) 
of the output SL equation \21} satisfies KdV equation ([l]); 



It 



3 , ,1 ,/ 



for all x,t > 0. 

Let us calculate the value of q{0,t),t > 0. Using formula (|4ip we obtain that 



g(0,t) = -2^[ 1 ] B''{x,t)X'\x,t)B{x) 



dx 



Using ([26} and ([T4)l it becomes 
-ig(0,t) =tr 



i 

-i 



B''{x,t)X-\x,t)B{x))\x=o - (-)'U=o 

T 



Since — = Jg q{y, t)dy at zero this term vanishes and we obtain 



--g(0,t)=tr 



i 

-i 



B*{x,t)X-\x,t)B{x,t))\x=o 
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Moreover, from (|39p and (|26p it follows that 



B'{x,i 



and as a result 
1 



q{Q,t) =(tr 
= tr([ 



i 
-i 
1 




,/>(0,t) 
r(0,f) 



B*(:E,i)X-i(a;,i)B(a;))U^o = tr 
:os((\/A)^t) -i\/Asm((\/A)^t) 



]* = [ 


1 " 














i 




1— ' 


— i 





[ 






= -/^^/Asm((^/A)3^)da;(A). 
in other words, it is assigned by the vessel construction: 



i 
-i 



qiO,t) 



VXsm{{V\ft)du(\). 



\*B{0,t) = 



(48) 



In order to satisfying different initial conditions for q{0, t) we have to perturb the original spectral 
function a;(A) as follows. Suppose that the measure i^(A) satisfies: 



3(\/Ax)cos(\/Ay)di/(A) = 0, Vx,?/>0. 



Let us define a new measure 



p'(A) = p(A) + !^(A), u' = oj~\-v. 



Then the formula ([8^ with p'(A) will produce the same /(a, y) and hence by the uniqueness of GL 
equation, the same K(x,y), and hence the same potential. But the formula for the initial value for 
5(0, t) will become 

q'{Q,t)^2 ( y/\sxa({^/Xft)duj'{\)^q{Q,t) + 2 ( y/\sm{{^/\ft)dv{\), 
Jm Jm 

creating a perturbation of the original potential. 

Finally we notice that the uniqueness of the solution of the KdV equation in this case "follows" 
from the fact that choosing a spectral function p(A) we actually uniquely determine a generalized 
vessel of the corresponding pre-vessel. As a result, the solution can be expected to be unique, as it 
is parallel to the solution of the KdV equation on the whole line. 

Finally, we obtain that evolving the operators as in the KdV case, we can create solutions of the 
KdV equation for generalized potentials using of pre-vessels. 

Definition 5.2. The KdV evolutionary pre-vessels is defined as follows 

pre2J = {A, B{x, t), X{x, t); ai, aa, 7; LI,C^), 
where A = ifi and B{x,t),X{x,t) are defined by (|46|l . (I47p and satisfy the regularity assumption 

(ESI. 

Theorem 5.3. Given a KdV pre-vessel, the generalized potential 7*(a;,t) exists for values x,t 
for which T{x,t) 7^ 0. At this point KdV equation ((T| holds. 
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6 Supplementary results 



6.1 Models of vessels in the case var[p{\)] < oo 

In the case var[p{X)\ < cxa it is known that the spectral function can be represented as follows 

oo 

dp{\) = f',{X)d\ + ^ bn5{X ~ A„), 

where /c(A) is an integrable function on R and 6„(5(A — A„) is the usual point measure, concentrated 
at A„ with the value fe„ there. We construct aj(A) from this function as at the Background 



pW, 



- VX, A > 0, 
A < 



that a{x) = 



s{VXx 



X 



-dw{X) possesses fourth continuous derivative. 



Define an element 
sequence {gn} as follows 

U = 

Define next 



5(A) 



G H to be a vector consisting of a function g{X) and an infinite 



A 



r. OO 

/ |3(A)|Vc(A)dA + ^|<7„|26„ <«D 

n 1=1 



5(A) 



iXg{X) 
{iX„g„} 



Bo 



{ 



} 



And apply the standard construction of a vessel to this data. Then 
B{x)-- 



[ cos{^/Xx) i\/Asin(\/A2;) ] 
{[ cos(vX^a;) sm{^X^x)} ] 



and 

X{x) = I+ r B{y)o2B\y)dy 
Jo 

will construct a vessel 93 using Theorem 13.41 Evolving of this vessel with respect to t will produce 



B{x,t) 



[ cos{^/Xx - {y/Xft) i^/Xsm{y/Xx ~ {y/Xft) ] 
{[ cos{^Kx{^%^ft) iy/Ksm{^%^x{^%^)^t)} ] 



and the corresponding X{x,t), defined by (|47|) producing a solution of the KdV equation ([T]) with 
the initial q{x) possessing the spectral function p(A) (var[p{X)] < oo). 

Notice that it can be generalized to a model involving (tempered) distributions, in which case 
the formula (fS^ will hold at the sense of distributions. This is a separate from this work project. 
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6.2 Construction of a vessel using Fadeyev theory 



The same approach can implemented for the setting of Fadeyev theory [Fad 74] . For this it is enough 
to consider Jost solutions (j){x,X) and tp{x,X) satisfying 



(f){s,x) — cos(sa::) + 
, , , sinfsa;) 



sin(s(a:: - t/)) 



<l{y)<t>{s,y)dy 



sm{s{x - y)) 



<i{v)i'{s,y)dy, 



for each s at the upper half plane, this actually means that these two functions satisfy ^ with the 
given q{x) and the initial conditions 

0(5, 0) = 1, (l)'{s, 0) = 0, i/'(s, 0) = 0, ip'{s, 0) = 1. 

Then define the following function, denoting I3{x) — q{y)dy 



S{\,x) 



(j){s,x) iip{s,x) 
-i{(f)'{s,x) — j3{x)(f){s,x)) ^'{s,x) — I3{x)tl}{s, x) 



cos{sx) isin(sx) 
ssin(s2;) cos(sa:) 



where we notice that all the functions actually depend on = Using the Fadeyev's work 

[Fad74) ■ or more precisely the work of Agranovich-Marchenko [MV60) it follows that this function 
is well defined for all x and satisfies the following differential equation 



^'(A, x) = a^\(Ji\ + 7,(a;))5(A, x) - ^(A, x)a^\ai\ + 7), 



^,{x) = 



for 7,(2;), defined by 

-i{P'{x)~l3\x)) 

I3{x) i 

At j: = it is a function of the complex variable A, which has jumps on the positive part of the 
imaginary axis and a finite number of isolated points on the negative imaginary part. In [Mela] it is 
shown that there exists a vessel with this transfer function with initial 5(A,0), where the operator 
A is again diagonal. Applying KdV vessel construction for this 5(A,0), we will obtain a solution of 
((ll, analogous to the classical solution using Fadeyev's inverse scattering. Notice that in this case 
the measure cj(A) satisfies var[ijj{\)] < 00 and we can use the standard construction to model the 
vessel. 



6.3 Examples of equations possessing vessel solutions 

The work in this section is presented in [Melb] for the case of bounded operators, but the ideas 
presented here are immediatly generalizable for an unbounded A, generating a Co semi-group. 

Definition 6.1. NLS vessel parameters are defined as follows: 



01 



And constructing a vessel in a similar to SL manner from a realized function 5(A,0), we will 
obtain solutions of the NLS equation: 





■ 


1 


■ 1 ■ 




■ 


■ 





1 




-1 


, 7 = 









^y{\,x) - {-iX)a2i\,x) + 



p{x) 
-p-ix) 



u{x, A) = 0, 
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where I3{x) is easily derived from the vessel. Evolving next the obtained vessel using (|43|l . (I44p . we 
will obtain solutions of the evolutionary NLS: 



For the Canonical systems, one can use the following parameters: 
Definition 6.2. Canonical systems vessel parameters are defined as follows: 



(49) 





■ i ' 




1 


■ 




■ 


■ 












1 


, 7 = 









and to use a similar theory of evolutionary vessels to create solutions of 



4/32 



where it turns out that for afixed if G R 7, (a;,t) 



'{x + Ky 



-2il3{x,t) ih{x,t) 
ih{x,t) 2il3{x,t) 



(50) 

with h{x, t) — 



1 



(x + K)^ ■ 
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